Transition between localized and extended states in the hierarchical Anderson model 



CO 

o 

CN 
u 

Or 

< 

CN 
CN 



CO 



c3 



i 

C 

O 
o 



CN 
> 
CN 

o 

CN 

o 



F. L. Mctz 1 , L. Lcuzzi 1 ' 2 , G. Parisi 1 ' 2 - 3 , V. Sacksteder IV 4 
1 Dip. Fisica, Universita La Sapienza, Piazzale A. Moro 2, 1-00185, Rome, Italy 
IPCF-CNR, DOS Roma Kerberos, Universita La Sapienza, P. le A. Moro 2, 1-00185, Rome, Italy 
3 INFN, Piazzale A. Moro 2, 00185, Rome, Italy 
4 Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 

(Dated: April 23, 2013) 

We present strong numerical evidence for the existence of a localization-delocalization transition 
in the eigenstates of the 1-D Anderson model with long-range hierarchical hopping. Hierarchical 
models are important because of the well-known mapping between their phases and those of models 
with short range hopping in higher dimensions, and also because the renormalization group can 
be applied exactly without the approximations that generally are required in other models. In 
the hierarchical Anderson model we find a finite critical disorder strength W c where the average 
inverse participation ratio goes to zero; at small disorder W < W c the model lies in a delocalized 
phase. This result is based on numerical calculation of the inverse participation ratio in the infinite 
volume limit using an exact renormalization group approach facilitated by the model's hierarchical 
structure. Our results are consistent with the presence of an Anderson transition in short-range 
models with D > 2 dimensions, which was predicted using renormalization group arguments. Our 
finding should stimulate interest in the hierarchical Anderson model as a simplified and tractable 
model of the Anderson localization transition which occurs in finite-dimensional systems with short- 
range hopping. 



I. INTRODUCTION 

After more than fifty years the Anderson transition^ 
between localized and extended wave-functions of a sin- 
gle quantum particle moving in a disordered medium 
remains the focus of considerable interest Crucial 
contributions to this field have been made by exactly 
solvable tight-binding models, such as 1-D models with 
nearest-neighbour hopping 2 ^ and models on the Bethe 
lattice££ Here we consider another interesting class of 
tight-binding models with long-range hopping arranged 
in a hierarchical block structure and decaying according 
to a power law with exponent a. Hierarchical models 
have a long history in statistical physics starting with 
Dyson^, and (as we will explain later) they provide an 
indirect route to understanding phases and critical be- 
haviour in Z)-dimensional systems^ 

Wc study the hierarchical Anderson model (HAM) in- 
troduced by Bovier, which combines on-site disorder with 
hierarchically-structured long range hopping^ In the ab- 
sence of disorder, the spectrum is an infinite set of highly 
degenerate flat bands that accumulate at the upper spec- 
tral edge. The degeneracies are arranged in a geometric 
series: one half of the pure HAM's states lie in the lowest 
energy band, one quarter in the next highest energy, etc. 
Hierarchical models preserve their structure under renor- 
malization group transformations r-~— which has allowed 
proof of several rigorous results about the site disordered 
HAM's spectrum;&i2r— and may promise exact exten- 
sions of the successful scaling theory of localization^ In 
particular, the absolutely continuous part of the spec- 
trum vanishes and the model presents only spectral lo- 
calization, provided that the hopping decays sufficiently 
quickly with distance, i.e., the hopping decay exponent 
a > 3/2^3. 



Unfortunately, much less is known about the size of the 
HAM's eigenvectors. The degeneracies of the pure model 
permit different choices of mutually orthogonal sets of 
eigenvectors. The most extended set consists of infinitely 
extended plane waves, while the least extended set has 
sizes that are strongly band-dependent, with very local- 
ized states in the lowest band and infinitely extended 
states in the highest band. In the presence of on-site 
disorder, it recently has been argued that all states are 
always localized^ based on an analogy with the criti- 
cality results for random-matrix models, such as ensem- 
bles of ultrametri o 18 ! 19 and power-law random banded 
matrices^ Both models, characterized by an exponent 
a controlling the power-law decaying random hoppings, 
exhibit an extended phase for a < 1 and a localized phase 
for a > 1. This would rule out the possibility of a tran- 
sition in the HAM, which has a well-defined macroscopic 
limit only for a > 1. However models with random hop- 
ping are relatively simple: the scattering length vanishes 
and only the localization length is important. The HAM 
belongs instead to the class of models with deterministic 
hopping, which are much richer because they have non- 
trivial physics at both length scales. In particular, the 
1-D Anderson model with on-site disorder and determin- 
istic power-law hopping exhibits a localization transition 
at its upper spectral edge^i - — 

In this work we show that the HAM exhibits a 
localization-delocalization transition near its upper spec- 
tral edge. We perform a thorough numerical study of 
the inverse participation ratio, which is the inverse of the 
eigenstate volume. Thanks to the HAM's invariance un- 
der block renormalization group (RG) transformations, 
we obtain recurrence equations for calculating the resol- 
vent matrix. This recursive method allows us to calculate 
the IPR in systems large enough to precisely determine 
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ratio P(E). The former is defined as 
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FIG. 1: Schematic representation of the hierarchical Anderson 
model, cf. Eq. JT]), with L — 2 3 sites. Lines denote hopping 
energies i p between sites in distinct blocks of size 2 P ~ . 



their infinite size behavior. Our results also suggest that 
there is a critical value of a above which all states arc 
localized, in analogy with the lower critical dimension 
D = 2 below which finite-dimensional short-range sys- 
tems are always localized and above which an Anderson 
transition was predicted using RG arguments*^ 



II. THE HIERARCHICAL ANDERSON MODEL. 

The HAM is a 1-D tight binding model with L = 2 N 
equally spaced sites, and independently distributed ran- 
dom site potentials ej, i = 1,...,L. The Hamiltonian 
reads 



n N = 5>|i><*| (l) 

N 2 N - p 1.2 p 

+ E Vp E E I ( r - !) 2P + *> ((»• - !) 2P + 3 I , 

where \i) is the canonical site basis. The second line is 
the hierarchical hopping matrix introduced by Dyson*£ 
It is the heart of the hierarchical Anderson model, and 
is organized in a tree as illustrated in Fig. [TJ The high- 
est level of the tree has index p = N and the lowest 
level has index p = 1. At each level the system is di- 
vided into 2 N ~ P separate blocks, each of which contains 
2 P sites. The hopping between any two sites within a sin- 
gle block has energy V p . As seen in Fig. [TJ the hopping 
between sites in two different blocks is determined by lev- 
els higher in the hierarchy and has energy t p = Y^n=p 
We study the deterministic HAM, which has hopping en- 
ergies V p = 2~ a ( p ~ 1 \ This exponential decay in the level 
index p ensures that in large N ^> 1 systems the hop- 
ping energy between sites separated by a distance O(L) 
decays according to a power law t p oc 0(L~ a ), the same 
as 1-D Anderson models with power-law hoppingj2i~— 

We study the infinite volume limit of the average den- 
sity of states (DOS) p(E) and of the inverse participation 
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p(E)= lim (-J2S(E-E I 

L—too \ Li z — * 



(2) 



where (...) is the average with respect to the disorder po- 
tential ti and E^ are the HAM's eigenvalues. The DOS 
measures the averaged spectrum, but does not contain 
any signal of the eigenstates' localization or derealiza- 
tion. We therefore study the average inverse participa- 
tion ratio (IPR) of the normalized eigenstates j?/^):— ~— 
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= lim 
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where l£ = J2i=i ((*lVv)) 4 is the IPR of an individual 
eigenstate. Its inverse measures the eigenstate's volume. 
The IPR is restricted to the interval < P(E) < 1. 
States that are perfectly localized on a single site satisfy 
P{E) = 1, and states that are equally distributed across 
all sites satisfy P(E) = l/L->0. 

In the pure W = HAM the DOS is a series of flat 
bands p puTe (E) = Y%Li 2- p S(E-E^). Each flat band 
is related to a level in the HAM's hierarchy. The bands' 
degeneracy decreases repeatedly by factors of two as one 
moves to higher energy, thus yielding the factor 2~ p . The 
difference between consecutive energetic levels falls off as 
E^-EP UTC oc 2"( Q " 1 )P and, hence, these accumulate at 
the upper spectral edge E££ Te . Near E^ lrc the integrated 
density of states Af(EV nie ) = Y7i=i follows a power 
law similar to that of short-range finite-dimensional sys- 
tems: A/-(£P urc ) = 1 - C (E^ c - EP ul - c ) ds/2 M^^. 
Here d s = 2/ (a — 1) is the spectral dimension which 
controls both diffusion and the long-distance physics 
of second-order phase transitions such as the Anderson 
transition. For W > 0, the integrated DOS of the HAM 
exhibits a Lifshitz tail at the upper spectral edge, with a 
Lifshitz exponent given by the spectral dimension^ This 
is the same behavior as observed in short-range finite- 
dimensional systems with on-site disorder, where the in- 
tegrated DOS exhibits a Lifshitz tail controlled by the 
Euclidean dimension^ Overall, as -E£^ re is approached, 
the spectral properties of the pure HAM become sim- 
ilar to short-range finite-dimensional systems. This is, 
thus, the most promising region for studying localization 
transitions. Much of the interest in hierarchical models 
originates in their mapping to short-range models whose 
Euclidean dimension is strictly related to the spectral di- 
mension (see, e.g., Ref. [3(| and Rcfs. therein). 



III. RENORMALIZATION EQUATIONS FOR 
THE RESOLVENT 

We obtain the DOS and IPR from the diagonal ele- 
ments of the resolvent matrix G^ N \z) = (z — Hn)" 1 , 
where z = E — in, and rj is a small positive rcgularizcr 
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that smooths our numerical results over an interval in 
the spectrum with width proportional to t]M& We use 
the following formulas iS&SLi&M 



P(E) 



1 

lim lim — V^/lmG?,- \z 
n ^o+ L^oo Lit ^— ' \ 



(4) 



P(E) = lim Mm — ^— £ (|Gf ^z)! 2 ) . (5) 

The HAM's hierarchical structure allowed us to develop 
a block RG approach which recursively calculates the re- 
solvent for one instance of the disorder. Our calculation 
has two phases: a sweep up the hierarchy, and then a 
sweep back down. At each step I of the sweep up we 
remove the basis states associated with one flat band, 
and calculate an energy- dependent effective Hamiltonian 
which acts in the reduced basis but exhibits the same 
poles found in the original full-basis Hamiltonian. This 
effective Hamiltonian retains the hierarchical form but 
its hopping energies {Vp^} and disorder potentials {/4^} 
are renormalised according to 
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(i-i) . (t-i) 
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p=l,...,N-£. 



N-l 



(6) 
(7) 



Hopping energies and disorder potentials at the begin- 
ning of the sweep up, V t 



CO) 



Tlip=i Vpi are those of the original hierarchical Hamilto- 
nian. After I — N steps we reach the top of the hierarchy 
and obtain a single site effective Hamiltonian with disor- 
der potential /z^. The resolvent of this Hamiltonian is 
simply G?0) = -1/^*°. We use this resolvent to begin 
the sweep back down, in which we progressively restore 
the original basis and recursively calculate the resolvent's 
diagonal elements in the restored basis: 
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FIG. 2: Comparison of the HAM's average DOS and IPR: 
numerical diagonalization (filled circles) vs. RG method (solid 
lines), with hopping decay exponent a = 7/4. The energies 
of the pure model's flat bands are marked with cyan vertical 
lines in the upper pane. The inset shows that when the IPR 
is small the RG method is sensitive to the spectral line width 
77; the dashed and solid lines were obtained with r\ — 0.01 and 
77 = 0.005 respectively. 
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with 7^ 1 = ^2i-i + M2i ±> ■ This procedure yields the 
diagonal elements of the resolvent in the original system, 
and its memory consumption and computational time 
grow only linearly with The derivation of Eqs. ([6][9]) 
is presented in App. [A] 



IV. RESULTS 

Fig. [2] compares the DOS and IPR calculated with 
our renormalization method (solid lines) and rj = 0.005 
to standard numerical diagonalization (filled circles) in 



a system of size L = 2 10 . The potential is generated 
from a Gaussian distribution with zero mean and stan- 
dard deviation W. Diagonalization results are averaged 
over 7V e = 10 3 disorder realizations and renormalisation 
results over Af e — 2 x 10 4 realizations. Fig. [2] shows 
excellent agreement between the two methods. 

The only important discrepancy is found in the IPR 
at small disorder W = 0.6, where the DOS falls precipi- 
tously. The observed discrepancy is explained by Fig. [2js 
inset, which compares results with two values of the regu- 
larization parameter: 77 = 0.01 and 0.005. The latter lies 
closer to the diagonalization results, which indicates that 
when the disorder is small the limit 77 — > + is reached 
only at 77 < 0.005. 

Fig. [2] also gives an overview of the DOS and IPR 
across the spectrum for a representative hopping decay 
exponent a = 7/4 and four different values of the dis- 
order strength W = 0.6, 1.0, 1.4, f.8. At small disor- 
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FIG. 3: Average DOS and IPR near the upper spectral edge 
in a large L — 2 23 system at several disorder strengths and 
a = 3/2, T) = 5 x 10" 4 and Me = 30. 



der W = 0.6 the average DOS is separated into several 
bands whose positions coincide with the pure system's 
flat bands. The associated minima in P(E) show that the 
eigenstates are bigger in the band centers and smaller at 
the band edges. When the disorder is increased the bands 
progressively blur together and P(E) steadily increases 
as the eigenstates become ever more localized. Fig. [3] 
shows the same behavior at a = 3/2 in systems of size 
L = 2 23 . Reaching such large sizes allows us to explore 
smaller ?y values and obtain detailed results about many 
bands near the upper spectral edge. Indeed, in order 
to obtain statistically significant results the spectral line 
width rj must considerably exceed the mean level spacing 
[Np{E)]-\ 

In general the IPR exhibits several local minima cor- 
responding to large states near the centers of the HAM's 
bands, and the global minimum lies near HAM's upper 
spectral edge. In order to verify the existence of extended 
states at finite W we focus on the asymptotic value of 
the global minimum of the IPR, P m i n (W), in the L^oo 
limit and for infinitesimal r\ — > + . The main graph 
in Fig. [4] summarizes our calculation of Pmin{W) for a 
particular hopping decay a = 3/2 and disorder strength 
W = 0.8 which lie close to the derealization transition. 
We display the IPR of a very large L = 2 27 system at 
three different values of 77. Statistical errors at smaller 77 
are larger because of 77's proximity to the level spacing. 
App. [B] includes a detailed discussion of these errors in 
the limit 77 — > + . In particular, we have checked that for 
L > 2 26 the IPR curves at fixed 77 do not change with L, 
which signals that they accurately represent the infinite 
volume limit. 

The IPR depends on 77, and as 77 — >• + the global 
minimum deepens and shifts toward higher energy. This 
effect is not significant at larger disorder W > 1.0, but at 




W=1.0 , E=5.480 
W=0.95 , E=5.465 
W=0.9 , E=5.445 
W=0.8 , E=5.405 




FIG. 4: The 77 — > + limit. The three curves show the average 
IPR at three values of 77 and a = 3/2, W = 0.8, L = 2 27 , and 
M e = 10. The global minimum decreases and shifts to higher 
energy. The inset shows the IPR versus 77 at the energy of 
the global minimum. Again a = 3/2, but now L = 2 28 and 
Me = 100. The solid lines are linear fits to P(rj) = P m i n + 677. 
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Wc 


7t2 
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X 2 /ndf 


ndf 


3/2 


0.684(7) 


2.68(7) 


0.080(2) 


0.96 
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7/4 


0.016(5) 


3.57(8) 


0.105(3) 


1.21 


5 



TABLE I: Values of the parameters and the \ 2 of the power 
law fit to the IPR data shown in Fig. [5] W c is the critical 
disorder where the derealization transition occurs, and 712 is 
a critical exponent. 



smaller disorder it forces us to use considerable care with 
the 77 — >• + extrapolation. The inset in Fig. @] displays 
our extrapolation to the limit 77 — > + at four weak dis- 
order strengths. At each W we find the energy E min (W) 
of the local minimum at the lowest value of the spectral 
width parameter employed, 77 = 10~ 5 , and then graph 
the IPR at that energy as a function of 77. Concerning 
uncertainty in E m i n (W), we have checked that it affects 
our results only slightly, and in any case can only cause 
an unduly careful overestimate of the IPR. The fitting 
curves in Fig. @]show that the IPR depends linearly on 77 
via P(rj) = P m in + br]. This allows us to determine very 
accurately the asymptotic global minimum of the IPR. 

The straight lines in Fig. [5js log-log plot are the cen- 
tral result of our work: strong numerical evidence that 
the minimum IPR P m i n (W) converges to zero according 
to a power law P m i n (W) = A(W — Wc) 71 ' 2 , similar to 
the power law observed in finite dimensional short-range 
systems At smaller disorder W < W c the HAM 

exhibits a delocalized phase. Table Q] reports the best 
fit parameters for two values of the hopping decay expo- 
nent a = 3/2,7/4. In both cases our data excludes the 
possibility that W c = 0. 
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FIG. 5: Power law behaviour of the minimum IPR P m in(W) 
near the critical disorder strength W c where it converges 
to zero. Solid lines represent the power-law fit P m in(W) = 
A(W — W c ) n ' 2 , with parameters from table [T] 



V. CONCLUSIONS. 

We have analyzed the DOS and the IPR of the hierar- 
chical Anderson model by means of a RG-based calcula- 
tion of the resolvent matrix, finding strong evidence for 
a localization-delocalization transition at finite disorder 
at a = 3/2 and 7/4. Since it has been proven rigor- 
ously that the absolutely continuous part of the spec- 
trum vanishes for a > 3/2^ our results indicate that 
spectral localization may not imply the existence of ex- 
ponentially localized eigenvectors. A study of the spatial 
decay of the resolvent elements should clarify this point 
and we expect our work to stimulate further research in 
this direction. Our results also indicate that the HAM 
differs from the 1-D tight-binding model with power-law 
hopping ; 23 ' 24 where all states are localized for a > 3/22^ 
Since the HAM's spectral dimension can be mapped to 
the spatial dimension of Anderson models with short- 
range hopping, we expect that an Anderson transition 
exists in the regime 1 < a < 2, with a ~ 2 playing a 
role analogous to the lower critical dimension. The pres- 
ence of extended states in one dimension is not exclusive 
to models with long-range hopping, but it has been also 
observed in systems with short-range hopping and cor- 
related on-site disorder Lastly we mention that our 
RG method can be used to compute off-diagonal elements 
of the resolvent, allowing determination of other relevant 
quantities such as the longest localization lengthj^ 
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Appendix A: Derivation of the renormalization 
equations 

We discuss here how to derive the RG equations (6-9) 
that are used in the main text. This calculation can be 
performed using only linear algebra, but we find it more 
convenient to use Gaussian integrals. We first rewrite 
{G k (z)}k=i,...,L, the diagonal elements of the resolvent 
G^ N '(z) = (z — Hn)" 1 , as the Gaussian integrals 



Gl N \z) = , 



J dtp (j>l cxp [C 



')] 



(Al) 



£ W (^l,...,2") 



where d<fi = J 



f^expfrW^,...,^)] 

2 N 

l -J2^ + wW (o, a * ; v, v.) 



3 = 1 



and [Lj 



T N V 



The 



function encodes the hierarchical hoppings: 



N 2 N -"I 




-1)2* +j 



We have introduced the simplified notation Xi t ...^ = 
X\,...,xj(. The function C^ N > {<j>i ... 2 N ~) nas the same 
form as the HAM's Hamiltonian and therefore preserves 
its formal structure under a RG transformation: a local 
term incorporating the random potential and a non-local 
hierarchical hopping term. We make a change of integra- 
tion variables 
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j-l±4>2j), j = l,...,2 



N-l 



which transforms the hierarchical term as follows: 



w (N) {<h,...,2»;Vi,..., Jf )=iV 1 



3=1 



(,•: a « : :i ; v :)• 



where V' = 2V p . This transformation allows us to explic- 
itly calculate the integrals over {V^ 1.7=1, ...,2 N - 1 in Eq. 
(|A1[) . halving the number of degrees of freedom. After 
performing the transformation and integration we obtain 
an equation which relates {G^ (2)}i=i,...,z for the orig- 
inal model with L sites to {G^ 1 (^)}t=i,...,L/2 f° r a 
model with L/2 sites, but with renormalizcd parameters. 
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Partitioning {G\ N \z)} i=1 2 N into two sectors (one 
for the even sites and another for the odd sites) , we obtain 
the following expressions 



= 2 Jd^ + dxp-e HiN - 1) ^ ± ) 

z/dV> + dV>-(V>+-^) 2 e g(JV " 1) ^ ± ) 
2 fd^+dip-e^-'H^) 



(A2) 
(A3) 



In the above expressions we have changed integration 



variables to dtp' 



niy — 1 _i_ 

dipj and we have defined 
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FIG. 6: Size effects on the IPR at a = 1.5, W = 0.9, 77 = 10" 
and A/" e = 10. 



+7 £ Cjvjvj 
3=1 



where the following quantities are complex valued: 
1 

°* = 2 
1 

2 
1 

2 



A, 



C, 



{p>2j-\ + V2j) + 2Vl, 



Since Eqs. (|A2[) and (| A3|) involve only simple Gaussian 
integrals with respect to ip~ 2 „_ 1 , these variables can be 
integrated out one by one. We map the resulting expres- 
sion to Eq. (|Aip for a system with 2 Ar ~ 1 sites, rcnormal- 
ized disorder p,^ 2JV _i which obeys Eq. (8) in the main 
text, and renormalised hopping potential V[ We 
obtain Eqs. (6-9) at the first RG step £ = 1 of the orig- 
inal model. Performing these steps recursively leads to 
the recurrence equations (6) and (7) shown in the main 
text: 
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Appendix B: Performing the 77 — > + limit 
numerically 

The rcgularization parameter 77 in the resolvent gives 
each eigenvalue a line width proportional to 77. This can 
be understood by analysing our equation for the DOS 



p(E) = lim lim — V ( 
77— s-o+ L— >oo Lit ~ — ' 



(JV) 



(Bl) 




FIG. 7: 77 dependence of the ratios r// Al,at, v (E) (open cir- 
cles) and a L ,M,ri(E)/pL,M,n(E) (filled circles). Al,at, v (E) 
is the approximate mean level spacing and <JL,Af,-q{E) is 
the standard deviation around plM,v(E)- Results were ob- 
tained using Eqs. (6-9) in the main text with N — 28 and 
jV c = 100. We set (W = 0.8, E = 5.405) for a = 3/2, and 
(W = 0.2, E = 3.532) for a = 7/4. These values were used to 
produce the left-most (smallest disorder) data point in Fig. 5 
of the main text. 



The right hand side of this equation is the limit 77 — > 
+ of a sum of Lorentzian functions with width 77 and 
centered at E. The Lorentzians quantify the distances of 
That's eigenvalues from the energy E. As 77 approaches 
the mean level spacing from above our observables will 
display larger and larger fluctuations, since our averages 
will include smaller and smaller numbers of eigenstates. 
If 77 is smaller than the level spacing then one obtains 
results which have no physical meaning. Accurate results 
for very small 77 are obtained only if the system size L and 
the number of samples Af e are large enough. Fig. |5]shows 
how this issue influences the IPR. We fix the number of 
samples Af e and spectral line width 77 and vary the system 
size. Convergence is obtained at L > 2 26 . 

If we define pL,M,r){E) as the average DOS of a finite 
though very large system, we can estimate the mean level 
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spacing Al^^ v (E) around E as 

We estimate the error at small rj by calculating 
crL,Ar, v (E)/pL,Ar, v (E) and rj/A L ^ t71 , where a Li j^ tV (E) is 
the standard deviation of pL,M,ni.E). Typical results are 
displayed in Fig. [71 When we decrease rj — > + we find 
monotonic growth in <tl,M.t){E) / Pl,a/.ti(E) and mono- 
tonic decay in 77/ Ax, ,jv>r For small enough 7? we reach a 
regime where a L ,M,v( E )/ Pl,U,v( E ) = 0(1), r]/A L .^ tn = 
0(1), and pL,j\f,r)(E) exhibits large fluctuations. We con- 



clude that the limit rj — > + is achieved, for practical 
purposes, when A^jy^ <C 77 <C 1, i.e. in very large sys- 
tems. Therefore we establish a sensible lower cutoff on 
77 by imposing a maximum value of the average DOS's 
relative error clM^E)/ pL,M,-n{E). 

The results for the minimum IPR displayed in Fig. 5 
of the main text were obtained by choosing 77 = 1CP 5 
as the lower cutoff when a = 3/2 and 77 e [10 - ', 10~ 5 ] 
when a = 7/4. This ensures that the relative er- 
ror <jl,N \n{E) I Pl,N \n{E) is restricted to the interval 
[10~ 2 , 10 _1 ], as can be seen in Fig. [7J 



1 P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

2 B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 
(1993). 

3 F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80, 1355 
(2008). 

4 E. N. Economou, Green's functions in quantum physics 
(Springer, Heidelberg, 2006). 

5 R. Abou-Chacra, P. W. Anderson, and D. J. Thouless, J. 
Phys. C: Solid St. Phys. 6, 1734 (1973). 

6 R. Abou-Chacra and D. J. Thouless, Journal of Physics C: 
Solid State Physics 7, 65 (1974). 

7 F. J. Dyson, Commun. Math. Phys. 12, 91 (1969). 

8 S. Molchanov, Proc. Lukacs Symposiu p. 179 (1996). 

9 A. Bovier, J. Stat. Phys. p. 745 (1990). 

10 G. Baker, Phys. Rev. B 5, 2622 (1972). 

11 Y. Meurice, J. Phys. A: Math. Theor. 40, R39 (2007). 

12 E. Kritchevski, Proc. Am. Math. Soc. 135, 1431 (2007). 

1 E. Kritchevski, "Hierarchical Anderson model" in "Proba- 
bility and mathematical physics: a volume in honor of S. 
Molchanov", vol. 42 (Amer. Math. Soc, 2007). 

14 E. Kritchevski, Ann. Henri Poincare 9, 685 (2008). 

15 S. Kuttruf and P. Moller, Ann. Henri Poincare 13, 525 
(2012). 

16 E. Abrahams, P. W. Anderson, D. C. Licciardello, and 
T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979), URL 
http : //link . aps . org/doi/10 . 1103/Phys RevLett .42 ■ 673 | 

17 C. Monthus and T. Garel, JSTAT p. P05005 (2011). 

18 Y. V. Fyodorov, A. Ossipov, and A. Rodriguez, JSTAT p. 
L12001 (2009). 

19 E. Bogomolny and O. Giraud, Phys. Rev. Lett. 106, 
044101 (2011). 

20 A. D. Mirlin, Y. V. Fyodorov, F.-M. Dittes, J. Quezada, 
and T. H. Seligman, Phys. Rev. E 54, 3221 (1996). 

21 C. Yeung and Y. Oono, Europhys. Lett. 4, 1061 (1987). 

22 A. Rodriguez, V. A. Malyshev, and F. Dominguez-Adame, 



J. Phys. A: Math. Gen. 33, L161 (2000). 

23 A. Rodriguez et ai, Phys. Rev. Lett. 90, 027404 (2003). 

24 A. V. Malyshev, V. A. Malyshev, and F. Dominguez- 
Adame, Phys. Rev. B 70, 172202 (2004). 

25 F. A. B. F. de Moura, A. V. Malyshev, M. L. Lyra, V. A. 
Malyshev, and F. Dominguez-Adame, Phys. Rev. B 71, 
174203 (2005). 

26 F. Wegner, Z. Physik B 36, 209 (1980). 

27 F. L. Metz, I. Neri, and D. Bolle, Phys. Rev. E 82, 031135 
(2010). 

28 F. Slanina, Eur. Phys. J. B 85, 361 (2012). 

29 W. Kirsch and F. Martinelli, Commun. Math. Phys. 89, 
27 (1983). 

30 M. Ibanez Berganza and L. Leuzzi, arXiv:1211.3991v2 
(2012). 

31 M. Janssen, Phys. Rep. 295, 1 (1998), ISSN 0370-1573. 

32 Y. Song, W. A. Atkinson, and R. Wortis, Phys. Rev. B 76, 
045105 (2007). 

33 E. N. Economou and M. H. Cohen, Phys. Rev. B 5, 2931 
(1972). 

34 Y. V. Fyodorov and A. D. Mirlin, Phys. Rev. Lett. 67, 
2049 (1991). 

35 J. Bauer, T.-M. Chang, and J. L. Skinner, Phys. Rev. B 
42, 8121 (1990). 

36 T.-M. Chang, J. Bauer, and J. L. Skinner, J. Chem. Phys. 
93, 8973 (1990). 

D. H. Dunlap, H.-L. Wu, and P. W. Phillips, Phys. Rev. 
Lett. 65, 88 (1990). 
38 F. C. Lavarda, M. C. dos Santos, D. S. Galvao, and 



:i!> 



B. Laks, Phys. Rev. Lett. 73, 1267 (1994). 



F. A. B. F. de Moura and M. L. Lyra, Phys. Rev. Lett. 
81, 3735 (1998). 
40 Mathematically Eqs. (JS| and ((7|) are equivalent to recursive 
computation of the Schur complement. 



